Fidelity of Quantum Teleportation through Noisy Channels 
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We investigate quantum teleportation through noisy quantum channels by solving analytically 
and numerically a master equation in the Lindblad form. We calculate the fidelity as a function 
of decoherence rates and angles of a state to be teleported. It is found that the average fidelity 
and the range of states to be accurately teleported depend on types of noise acting on quantum 
channels. If the quantum channels is subject to isotropic noise, the average fidelity decays to 1/2, 
which is smaller than the best possible value 2/3 obtained only by the classical communication. On 
the other hand, if the noisy quantum channel is modeled by a single Lindblad operator, the average 
fidelity is always greater than 2/3. 
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Quantum teleportation Jl], [^] is a process by which 
a sender, called Alice, transmits an unknown quantum 
state to a remote recipient, called Bob, via dual clas- 
sical and quantum channels. Here a pair of maximally 
entangled particles, forming a quantum channel, should 
be used for the perfect quantum teleportation. However, 
while being distributed and kept by Alice and Bob, an en- 
tangled state may lose its coherence and become a mixed 
state due to the interaction with its environment. 

Bennett et al. Q noted that the quantum channel 
which is less entangled reduces the fidelity of telepor- 
tation, and/or the range of states that can be accurately 
teleported. Popescu || investigated the relations among 
teleportation, Bell's inequalities, and nonlocality. ft was 
demonstrated that there are mixed states which do not 
violate any Bell type inequality, but still can be used 
for teleportation. Horodecki et al. showed that any 
mixed two spin-i state which violates the Bcll-CHSH in- 
equality is useful for teleportation. Also Horodecki et 
al. Q proved the relation between the optimal fidelity 
of teleportation and the maximal singlet fraction of the 
quantum channel. Banaszek investigated the fidelity 
of quantum teleportation using non-maximally entangled 
states. Ishizaka [|| studied the quantum channel sub- 
ject to local interaction with two- level environment. Al- 
though the studies cited above reveal the important rela- 
tions between the degree of entanglement of the quantum 
channel and quantum teleportation, there seem to be lit- 
tle studies on the direct connection between the quantum 
teleportation and decoherence rates. Thus it might be 
interesting to know how the type and strength of noise 
acting on quantum channels affect the fidelity of quantum 
teleportation. 

In this paper, we investigate quantum teleportation 
through noisy channels by solving analytically and nu- 
merically a master equation in the Lindblad form. We 
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obtain the fidelity of quantum teleportation as a function 
of decoherence time and angles of an unknown state to 
be teleported. Thus we explicitly demonstrate Bennett 
et aVs argument that noisy quantum channels reduce 
the range of states to be accurately teleported. We also 
examine the characteristic dependence of the average fi- 
delity on types of noise acting on qubits at each stage of 
the teleportation. 

Let us consider quantum teleportation through noisy 
channels as illustrated in Fig. [j]. The top two qubits are 
taken by Alice and the bottom qubit is kept by Bob. Here 
measurements are performed at the end of the circuit for 
computational convenience. Classical conditional opera- 
tions can be replaced with corresponding quantum con- 
ditional operations || . Decoherence of an open quantum 
system is due to the interaction with its environment. 
Under the assumption of Markov and Born approxima- 
tions and after tracing out the environment degrees of 
freedom, the dynamics of an open quantum system is 
described by a master equation for the density opera- 
tor of the quantum system alone, p{t) in the Lindblad 

form Ga in 
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the i — th qubit and describes decoherence. Throughout 
this paper, <t£ denotes the Pauli spin matrix of the i-th 
qubit with a = x,y, z. The decoherence time is approx- 
imately given by 1/k^q. By switching on and off Ki i0l (t) 
we could control noise. We take the Hamiltonian of a 
qubit system as an ideal model of a quantum computer 
which is given by |12| 
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solid state qubits, various types of the coupling between 
qubit i and j are possible such as the XY coupling given 
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above, the Heisenberg coupling, and the Ising coupling 
JijOz Oz m NMR. The various quantum gates in Fig. [j] 
could be implemented by a sequence of pulses, i.e., by 
turning on and off B^(i) and Jij(t)- We develop the 
simulation code which solves Eq. (pi), the set of differ- 
ential equations for the density matrix p mn (t), based on 
the Runge-Kutta method Eq. ([[]) shows Tr p(t) = 1 
at all times. 

An unknown state to be teleported can be written as 
IV'in) = Qf|0) + /3|1) with \a\ 2 + \(3\ 2 = 1. It is convenient 
to rewrite |i/>i n ) Bloch vector on a Bloch sphere 

IVO =cosQ)e^ 2 |0)+sin^ e -^/ 2 |l), (3) 

where 9 and <p are the polar and azimuthal angles, re- 
spectively. The maximally entangled state of two spin-^ 
particles shared and kept by Alice and Bob is given by 

|/3oo) = ^(|00) + |11)). (4) 

The input state of the quantum teleportation circuit in 
Fig. [I] is the product state of \ipi n ) and |/3oo)- After the 
implementation of the quantum circuit of Fig. [l] and the 
measurement of the top two qubits, Bob gets the tele- 
ported state |Vw)- It is useful to describe the telepor- 
tation in terms of density operators 

Pout = Tri )2 {J7telPin ® Pen^tel}' ( 5 ) 

where p in = \ip ia ) (ip in \, p en = |/3 00 ) (f3 00 \, and Tri >2 is a 
partial trace over qubits 1 and 2. The unitary operator 
Utei is implemented by the teleportation circuit as shown 
in Fig. [lj. If the teleportation is ideal, the density matrix 
tcleportcd p on t is identical to pi n up to the normalization 
factor. 

As illustrated as dotted boxes in Fig. [j], we consider 
four different noisy channels, A, B, C, and D. In case A 
an unknown state \ipi a ) loses its coherence and becomes 
a mixed state before it is teleported. In case B an en- 
tangled pair, forming a quantum channel, become noisy 
while being shared and kept by Alice and Bob. In cases 
C and D, while Alice and Bob perform the Bell measure- 
ment and the unitary operation, respectively, noise may 
set in. For cases A and B we obtain both analytic and 
numerical solutions of Eq. ([!]). While in cases C and D 
the numerical solutions of Eq. (Q) are obtained. For our 
numerical calculation, Ki i0 ,(t) is turned on for the time in- 
terval r corresponding to the width of each dotted boxes 
in Fig. 0. 

The properties of quantum teleportation through noisy 
quantum channels are quantified by the fidelity which 
measures the overlap between a state to be teleported 
IV'in) and the density operator for a teleported state p out 

F(6,4>) = (il) in \p oat \ip in ). (6) 

Here the fidelity F(9, 4>) depends on an input state as 
well as the type of noise acting on qubits. We calculate 



F(6, 4>) and determine the range of states |^i n ) which can 
be accurately teleported. Since in general a state to be 
teleported is unknown, it is more useful to calculate the 
average fidelity given by 

Fw = —j dO J deb F '(0,4,) sin 9, (7) 

where Ait is the solid angle. 

Case A: states to be teleported are mixed. Alice is not 
able to know or copy the state to be teleported without 
disturbing it. So it may be pure or mixed. As Bennett et 
al. noted, the linear property of quantum teleporta- 
tion enables one to teleport not only a pure state but also 
a mixed state. The quantum operation £ transforms a 
pure state p m = \ipia) (ipm\ to a mixed state £ (pin)- The 
time-evolution of pure states to mixed states is described 
by Eq. (|IJ). See Ref. p|for the connection between two ap- 
proaches. From Eq. (pf), quantum teleportation of mixed 
states reads 

£(A>ut) = ^ ia {U tcX £(p- m )® Pcn ul x ) . (8) 

The decoherence of the state to be teleported, £(pi a ) is 
transfcred to the state teleported, £(p ut)- For various 
types of noise, we obtain both analytic and numerical 
solutions of Eq. ([!]) , and calculate the fidelity. 

Suppose a state to be teleport is subject to the noise 
Li iZ . It is easy to find the analytic solution of Eq. (|l|) 
when Hs(t) = 0. We obtain the mixed state to be tele- 
ported, £(p in ) as p<°°>(t) = Pl ( ° 0) (0), p^(t) = p^(0), 
and p( 01 )(t) = p i ( ° 1) (0)exp(-2Kt). Then from Eqs. (|) 
and (0), the fidelity can be calculated as 

F(M) = l-i(l-e- 2Kr )sin 2 #. (9) 

If 2kt < 1, F(6, 4>) ^ 1 - kt sin 2 9. On the other hand, 
if 2kt > 1, F(9, 4>) ~ ±(1 + cos 2 6). Fig. | (a) is the plot 
of Eq. H for 2kt = 3.0. 

Let us consider the state \ipi n ) is subject to the noise 
described by L\. x . After some calculations, we obtain the 
fidelity 

F(9,(j)) = i[l + sin 2 6»cos 2 (/) 

+ e- 2KT (cos 2 6> + sin 2 6>sin 2 0)] . (10) 

If 2kt < 1, F(9, cj)) ~ 1-kt(cos 2 6»+sin 2 9 sin 2 cp) . In the 

limit of 2kt > 1, we have F(9, 4>) ~ \(l + sin 2 9 cos 2 (f>). 
The plot of Eq. @ at 2kt = 3.0 is shown in Fig. | (b). 

Substituting Eqs. (^Tj) or (|) into Eq. (^), we get the 
average fidelity 

^v(r) = ^ + i e - 2 ^. (11) 

In Fig. |^, the solid line (denoted by Case A-l) shows the 
plot of Eq. (0), the average fidelity as a function of kt 
for the noise modeled by L± x or L\ z . 
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Now suppose the isotropic noise (L± x ,Lxy, and L% z ) 
is applied to the state |"0in) - The analytic solution of 
Eq. (jy) gives us the fidelity written by 



F( 



1 1 

= 2 + 2 e 



4kt 



(12) 



If 4kt < 1, F(6,<f>) ~ 1 — 2kt. For Akt > 1, we have 
F(8, <f)) ~ 5 as shown in Fig. | (c). In Fig. |, the dotted 
line (denoted by Case A-2) is the plot of Eq. (|l2|). 

Case B: Quantum channels are noisy. While being 
distributed and stored by Alice and Bob, an entangled 
state of two spin-i particles may be subject to noise. The 
dynamics of an entangled pair subject to quantum noise 
is described by the quantum operation £ acting on the 
pure entangled state, p on — > £(p cn ) or by Eq. (|lj). From 
Eq. (|^), the quantum teleportation with noisy quantum 
channels can be written as 

£(/w) - Tri, 2 {C/ te] Pin <g> £( Pon ) U f tcl } . (13) 

Wc find that the quantum teleportation process transfers 
the decoherence of the entangled pair £ (p en ) to that of 
the output state £(p ou t)- It should be noted that the 
quantum operation acting on the entangled pair £(p cn ) 
is a 4 x 4 matrix but effectively a 2 x 2 matrix. Thus 
overall features of case B are similar to case A except 
decoherence rates. 

Consider the quantum channel subject to the noise 
acting in one direction, for example, the z direction. 
This type of noise is modeled by Lindblad operators, 
L2, z = y/K2,z Oz and L 3 , z = ^k 3 . z o-f\ acting on 
an entangled pair, qubit 2 and qubit 3, respectively. 
Here we assume the same strength of decoherence rates, 
k = K2, z = Ks tZ . We obtain the fidelity F(9,4>) with the 
same form of Eq. ([)]) except the replacement of 2kt with 
4kt. That is F(9, ft) = 1 - |[1 - exp(-4Kr)] sin 2 9. For 
the noise described by L^ tX and L 3 x , the fidelity F(9, ft) 
is identical to the form of Eq. ([ll]) with exponent 4kt. 

Let us discuss Bennett et aZ.'s argument: the imperfect 
quantum channel reduces the range of state IV^m) that is 
accurately teleportcd 0|. Fig. || shows the fidelity F(9, ft) 
for various types of noise at 4kt = 3.0. For the noisy 
channel defined by L 2 , z and L 3 ^ z , the fidelity F(9 1 ft) is 
always the maximum value 1 at 9 = 0, 7r irrespective of 
kt as depicted in Fig. || (a). These angles indicate states 
|0) and |1), which are eigenstates of <r z . From F(9, 0) ~ 
|(1 + cos 2 9) in the limit of 4kt ^S> 1, the range of states 
to be teleported with fidelity F > 3/4 is determined by 
< 6 < tt/4 and 3?r/4 < 9 < n. The teleported states 
with fidelity 2/3 are in the region determined by cos 9 > 
l/y/3 or cos 9 < When Li, x and L 3 ^ x are applied 

to the qubits 2 and 3, we get F{9,ft) = 1 at 6 = ir/2 
and (f> = 0,tt for 4kt 3> 1, which shown in Fig. § (b). 
These angles represent states = ^(|0) + |1}) and 

= t=(|0) — |1)), i.e., eigenstates of a x . The range 
of states accurately teleported is depicted by contours in 
Fig. | (a) and (b). 



When the quantum channel is subject to noise in an 
one direction, we obtain the average fidelity as depicted 
in Fig. | (denoted by Case B-l) 



F av (r) 



-4kt 



(14) 



The average fidelity decays exponentially to the limiting 
value of 2/3. This is the best possible score when Al- 
ice and Bob communicate each other only through the 
classical channel ||] . 

Consider the case that the quantum channel is affected 
by isotropic noise, which is described by six Lindblad 
operators, £ 2 .q, and L 3 a with a = x,y, z. Then the 
analytic calculation of the fidelity can be written by 



F w = = - + 



1 1 

2 + r 



(15) 



As depicted in Fig. || (c), the fidelity F(6, 4>) is indepen- 
dent of angles of input states, 9 and <p f° r any value of 
kt. For the quantum channel subject to isotropic noise, 
one could not find the range of states that is accurately 
teleported. As shown in Fig. || (Case B-2), the average 
fidelity decay to the value 1/2. The number 1/2 can be 
obtained when Alice and Bob can not communicate at 
all and Bob merely selects a state at random. 

It should be noted that except decoherence rates kt, 
the overall features of cases A and B are identical. This 
implies that if a state to be teleported is realized by a 
single particle not an ensemble, one may not be able to 
identify whether the state to be teleported is mixed or 
the quantum channel is noisy. 

Cases C and D: Noise during Bell's measurement or 
the unitary operation. When Alice performs the Bell's 
measurement or Bob does the unitary operation on his 
particle of an entangled pair, noise may take place as 
depicted by the boxes C or D in Fig. [[} In contrast 
to cases A and B, it seems to be difficult to find an- 
alytic solutions of Eq. (|l]) for cases C and D because 
of the time-dependence of the qubit Hamiltonian Hs(t). 
Alice's Bell measurement on qubits 1 and 2 could be 
done by a controlled-not gate (CNOT) on qubits 1 and 
2, and a Hadamard gate Hi on qubit 1 as shown in 



Fig. 
Eq. (@ 



With a qubit system modeled by Hamiltonian 
the CNOT gate acting on qubits 1 and 2 could 



be implemented by the pulse sequence 



e-^R^i^R^i-^U^mRUn^mil,. Here 



nit 



CNOT = 



R 



. ]x {9) = e I<Ta = ' is a rotation of qubit j by angle 9 
about the x axis. A two qubit operation U^(9) on qubits 
1 and 2 is implemented by turning on the coupling J12 
for a time t corresponding to 9 = Ji2t/h. During each 
qubit operation, the noise modeled by Lindblad opera- 
tors is also switched on. Thus it seems to be not simple 
to obtain an analytic solution and we takes a numerical 
method to solve the problem. 

Consider the noise modeled by the Lindblad opera- 
tors, Li z and L 2z for case C, and L 3z for case D. Here 
the noise is switched on during the time interval r cor- 
responding to the total operation time which it takes to 
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implement Bell's measurement or controlled X and Z op- 
erations. The time interval r depends on the operation 
times of a single gate or a two qubit gate, proportional 
to /i/|B^)| and h/Jij, respectively. Fig. |I] shows the fi- 
delity F(9, </>) as a function of angle 9 for various values 
of kt. In contrast to the previous cases (Case A-l, Case 
A-2, Case B-b and Case B-2) whose fidelity is given by 
Eqs. (H) or ([h]), in cases C and D the degrees of the de- 
pendence of fidelity F(6, 4>) on angles 9 is maximum at a 
certain value of kt. Fig. [| (c) shows the differences be- 
tween the maximum and minimum values of the fidelity, 
q(kt) = m&x{F(9, </>)} — mm{F(9,<fr)}. It is not clear 
why <?(kt) has the maximum at kt « 0.98. As depicted 
in Fig. ||, the average fidelity falls to the value 1/2 and is 
approximately fitted by 

FAr) = \ + \e~ 1 - 2bKT ■ (16) 

One sees that cases A and B entirely differ from cases 
C and D. Although an analytic solution for case C and 
D can not be obtained, it can be understood why the 
average fidelity decays to 1/2 despite noise described by 
the Lindblad operator acting in one direction. Consider 
a rotation of a qubit about the x axis in the presence of 
noise modeled by L z . A simple calculation shows that 
the Bloch vector r of a qubit. which is defined by p — 
i(l + r-cr), falls to zero for any initial state. This means 
the qubit is depolarized and becomes a totally mixed 
state. Thus the average fidelity decays to 1/2 when a 
gate operation is done in the presence of noise. 

It is valuable to discuss our results in connection with 
the previous studies || . In Ref . ||, Popescu illustrated an 
example of a mixed pair which does not violate any Bell 
inequality but has an average fidelity 3/4 for arbitrary 
input states, given by p = |/ + \ (*~| with |*~) = 



-^|(|01) — 1 1 0) ) . Our calculation in which the quantum 
channel is described by this mixed state shows the fidelity 
F{6,4>) = 3/4, independent of angles 9 and <f> and thus 
gives us an average fidelity F av = 3/4. Horodecki et. 
al. H showed that the optimal fidelity of the standard 
quantum teleportation is given by / = (2Fab + l)/3, 
where Fab is the singlet fraction of the quantum channel. 
From Eq. (|l|), one can write Fab = (1 + e~ 4KT )/2. 

In conclusion, we calculated the fidelity and the aver- 
age fidelity of quantum teleportation subject to various 
types of noise during different steps of the teleportation. 
We examined the range of states that can be accurately 
telcported. Among states to be teleported, the eigenstatc 
of the Lindblad operators is less sensitive to the noise. It 
was shown that one can not distinguish whether an un- 
known state to be teleported, which is realized by a single 
particle, is mixed or the quantum channel is noisy. We 
found the dependence of the average fidelity on the type 
of noise affecting the quantum channel. If the quantum 
channel is subject to isotropic noise, the average fidelity 
may decay to 1/2. On the other hand, if the noisy quan- 
tum channel is described by a single Lindblad operator, 
the average fidelity is always greater than the value 2/3, 
the best possible value which can be obtained only by 
the classical communication. 
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FIG. 1: A circuit for quantum teleportation through noisy 
channels. The two top lines belong to Alice, while the bottom 
one to Bob. M represents measurement. The dotted boxes, 
A, B, C, and D denote noisy channels. Time advances from 
left to right. During the time interval corresponding to the 
width of the dotted box, the Lindblad operator is turned on. 




FIG. 2: Fidelity F(0, <j>) as a function of angles 8 and <j> of 
the state to be teleported for case A at 2kt — 3.0 and for 
case B at 4kt = 3.0. For case B in (a) the Lindblad operators 
Z/2,z and -L3, z, and in (b) L2,x and L^ tX are turned on. The 
maximum values of the fidelity of (a) and (b) are 1 and the 
minimum 1/2. The contours on the — <f> planes in (a) and 
(b) join the points with the fidelity F(0,<f>) = 3/4 and 2/3, 
respectively. In (c) the isotropic noise is applied. 
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FIG. 3: Average fidelity F^ v as a function of kt for cases A 
and B. The solid line (Case A-l) is the plot of Eq. (JTTJ) for the 
noise described by L\ x (or by Li, z ). The dotted line (Case 
A-2) is based on Eq. (|l|) corresponding to the isotropic noise. 
The dashed line (Case B-l) is for Eq. (|l4|), the noise modeled 
by 1/2, x and L^ :X (or by L2, z and £3,2). The dash-dotted line 
(Case B-2) is the plot of Eq. ( |l5| ) for the isotropic noise. The 
horizontal dotted line with 2/3 shows the maximum fidelity 
obtained only by the classical communication. 
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FIG. 4: Fidelity F(8,<j)) vs angle 9 for various values of kt. 
(a) for case C noise is modeled by L\ z and Z/2 Z , and (b) for 
case D by L^ z . F(8,<j>) is independent of angle <j> because of 
the cylindrical symmetry of Li Z with i = 1,2,3. Differences 
between the maximum and minimum values of the fidelity 
F(8, cj>) axe plotted as a function of kt. 



